Computer-aided simulation method for atomic-resolution scanning seebeck microscope (SSM) images

ABSTRACT

A computer-aided simulation method for an atomic-resolution scanning Seebeck microscope (SSM) image is provided. In the computer-aided simulation method, a computer may calculate a local thermoelectric voltage for a position of a voltage probe, to acquire an SSM image corresponding to the position.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a Divisional of co-pending U.S. patent application Ser. No. 14/299,170, filed on Jun. 9, 2014, which claims the benefit of Korean Patent Application No. 10-2013-0066020 and of Korean Patent Application No. 10-2014-0055617, respectively filed on Jun. 10, 2013 and May 9, 2014, in the Korean Intellectual Property Office, the disclosures of which are hereby incorporated herein by reference in their entireties.

FIELD OF THE INVENTION

The present invention relates to a computer-aided simulation method for an atomic-resolution scanning thermoelectric microscope.

BACKGROUND OF THE INVENTION

To observe a surface of a material in an atomic resolution, a scanning tunneling microscope (STM) technology disclosed in U.S. Pat. No. 4,343,993 may be used. The STM technology may be used to display an atomic-resolution image representing a surface of a material, based on a vacuum tunneling current localized by applying an external voltage to a gap between the material and a probe. However, in the STM technology, since the external voltage needs to be applied, Fermi electrons may be likely to be disturbed by the external voltage, and a narrow region may be scanned. Additionally, a sharp probe is required to generate a localized current, however, a yield may not be high due to a difficulty of forming a sharp probe.

To avoid the above disadvantages of the STM technology, a scanning thermoelectric microscope may be used instead of an STM. However, since it is generally known that it is difficult to localize heat in space, it may be difficult to actually implement the scanning thermoelectric microscope.

Recently, a technology of observing a surface of a material at an atomic resolution using a scanning thermoelectric microscope has been disclosed in “Thermoelectric imaging of structural disorder in epitaxial graphene” prepared by S. Cho, et al. and published in 2013 in arXiv: 1305.2845 (http://arxiv.org/abs/1305.2845) and Nature Mater. 12, 913.

Accordingly, it is desired to acquire an atomic-resolution image from a surface of a material, using a scanning thermoelectric microscope that is an actual apparatus, and to simulate a form of the surface using a computer to determine an atomic structure corresponding to the acquired atomic-resolution image.

SUMMARY OF THE INVENTION

An aspect of the present invention provides a computer-aided simulation method for an atomic-resolution scanning Seebeck microscope (SSM) image, to provide a method of simulating a form of a surface of a material corresponding to an image acquired by an actual scanning thermoelectric microscope and of recognizing the form.

According to an aspect of the present invention, there is provided a computer-aided simulation method for an atomic-resolution SSM image, wherein a computer calculates a local thermoelectric voltage V(r) for a position of a voltage probe, to acquire an SSM image corresponding to the position, using the following equation:

${V(r)} = {V_{diff} + {{S_{coh}(r)}{\int{{{\nabla{T\left( {r;r^{\prime}} \right)}} \cdot \frac{r^{\prime} - r}{{{r^{\prime} - r}}^{3}}}{\mathbb{d}^{3}r^{\prime}}}}}}$

In the equation, V_(diff) denotes a thermoelectric voltage drop in a diffusive transport region in a tip and a sample, S_(coh)(r) denotes a position-dependent Seebeck coefficient, r denotes a distance measured from a point voltage probe, r′ denotes material internal coordinates, ∇T(r; r′) denotes a temperature gradient radially weighted by a factor of 1/r², and

$\int{{{\nabla{T\left( {r;r^{\prime}} \right)}} \cdot \frac{r^{\prime} - r}{{{r^{\prime} - r}}^{3}}}{\mathbb{d}^{3}r^{\prime}}}$ denotes a volume integral of a temperature profile.

The volume integral

$\int{{{\nabla{T\left( {r;r^{\prime}} \right)}} \cdot \frac{r^{\prime} - r}{{{r^{\prime} - r}}^{3}}}{\mathbb{d}^{3}r^{\prime}}}$ may be defined by the following equation:

${\int{{{\nabla{T\left( {r;r^{\prime}} \right)}} \cdot \frac{r^{\prime} - r}{{{r^{\prime} - r}}^{3}}}{\mathbb{d}^{3}r^{\prime}}}} - {\Delta\;{T_{coh}(r)}}$

In the equation, ΔT_(coh)(r) denotes an effective temperature drop at an interface between the tip and the sample.

The position-dependent Seebeck coefficient S_(coh)(r) may be calculated from the Landauer formula by the following equation:

${S_{coh}(r)} = {{- \frac{1}{e\; T}}\frac{\int{{\tau^{el}\left( {E,r} \right)}\left( {E - E_{F}} \right)\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}{\int{{\tau^{el}\left( {E,r} \right)}\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}}$

In the equation, e denotes an electron charge, T denotes an absolute temperature, τ^(el) denotes an electron transmission probability, E_(F) denotes a Fermi energy, and ƒ denotes a Fermi-Dirac distribution function at the absolute temperature T.

The position-dependent Seebeck coefficient S_(coh)(r) may be derived to the following equation: S _(coh)(r)=S ^(tip) +S ^(sample)(r)

In the equation, S^(tip) denotes a tip Seebeck coefficient, and S^(sample)(r) denotes a sample Seebeck coefficient of a material.

The sample Seebeck coefficient S^(sample)(r) may be defined by the following equation:

${S^{sample}(r)} = {{- \frac{1}{e\; T}}\frac{\int{{N_{e}^{sample}\left( {E,r} \right)}\left( {E - E_{F}} \right)\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}{\int{{N_{e}^{sample}\left( {E,r} \right)}\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}}$

In the equation, e denotes an electron charge, T denotes an absolute temperature, E_(F) denotes a Fermi energy, ƒ denotes a Fermi-Dirac distribution function at the absolute temperature T, and N_(e) ^(sample)(E,r) denotes a local density of states (DOS) obtained from first-principles calculations for a material surface.

The effective temperature drop ΔT_(coh)(r) may be calculated by the Following equation:

${\Delta\;{T_{coh}(r)}} = {\frac{\Delta\; Q\; 3h}{\pi^{2}k_{B}^{2}T}\frac{1}{\tau^{ph}(r)}}$

In the equation, k_(B) denotes a Boltzmann constant, and τ^(ph)(r) denotes a phonon transmission probability of a tip-sample junction.

The effective temperature drop ΔT_(coh)(r) may be deduced by the following equation:

${\Delta\;{T_{coh}(r)}} = \frac{{V^{\exp}(r)} - V_{diff}}{S_{coh}^{theory}(r)}$

In the equation, Γ^(expt)(r) denotes an experimental thermoelectric voltage measured in a defect-free region, and S_(coh) ^(theory)(r) denotes a theoretical Seebeck coefficient of a material.

The effective temperature drop ΔT_(coh)(r) may be deduced by a function represented by the following equation: ΔT _(coh)(r)=ƒ(E ^(vdw)(r))

In the equation, ƒ(E^(vdw)(r)) denotes a function including a linear fitting formula or an exponential function, and E^(vdw)(r) denotes a van der Waals (vdW) energy The function may be based on a correlation between the effective temperature drop ΔT_(coh)(r) and the vdW energy E^(vdw)(r).

Effect

According to embodiments of the present invention, a computer may simulate an atomic-resolution image of a scanning Seebeck microscope (SSM) corresponding to a surface of a material based on a Seebeck effect.

BRIEF DESCRIPTION OF THE DRAWINGS

These and/or other aspects, features, and advantages of the invention will become apparent and more readily appreciated from the following description of exemplary embodiments, taken in conjunction with the accompanying drawings of which:

FIG. 1 schematically illustrates an atomic resolution of a scanning thermoelectric microscope according to an embodiment;

FIGS. 2A through 2F illustrate an experimental thermoelectric voltage image and a Seebeck coefficient image for a defect-free region of n-doped free-standing graphene, line profiles, and a correlation between an effective temperature drop and van der Waals (vdW) energy according to an embodiment;

FIGS. 3A through 3D illustrate images of a thermoelectric voltage for a point defect, a topographic image, and profiles according to an embodiment;

FIGS. 4A through 4F illustrate simulated images of Seebeck coefficients and thermoelectric voltages near point defects in n-doped free-standing graphene according to an embodiment;

FIGS. 5A through 5D illustrate a model and a density of states (DOS) of pristine free-standing graphene, vdW energy, and a vdW topography according to an embodiment;

FIGS. 6A through 6F illustrate locally-averaged Seebeck coefficients, a correlation between vdW energy and an effective temperature drop, and thermoelectric voltages according to an embodiment;

FIGS. 7A through 7D illustrate a model and a DOS of defective free-standing graphene, vdW energy, and a vdW topography according to an embodiment;

FIGS. 8A through 8D illustrate a model and a DOS of defective free-standing graphene with a single substitutional nitrogen (N_(c)), vdW energy, and a vdW topography according to an embodiment;

FIGS. 9A through 9D illustrate a model and a DOS of defective free-standing graphene with a defect complex (V_(c)-O_(c)) of a carbon vacancy (V_(c)) and substitutional oxygen (O_(c)), vdW energy, and a vdW topography according to an embodiment; and

FIGS. 10A through 10D illustrate a simulated thermoelectric voltage image for a V_(c)-O_(c), a fast Fourier transform (FFT) image corresponding to the simulated thermoelectric voltage image, an experimental thermoelectric image, and an FFT image corresponding to the experimental thermoelectric image according to an embodiment.

DETAILED DESCRIPTION OF THE INVENTION

Reference will now be made in detail to exemplary embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the like elements throughout. Exemplary embodiments are described below to explain the present invention by referring to the figures.

Hereinafter, in the present disclosure, a scanning Seebeck microscope (SSM) may refer to a scanning thermoelectric microscope based on a Seebeck effect. According to an embodiment, a computer-aided simulation may be performed on a material surface, to search for a material corresponding to an image of the material surface acquired by the SSM.

A scanning tunneling microscope (STM) may use tunneling, and may be used to acquire a direct wave function image based on a charge-transfer gap and a voltage bias. Unlike the STM, an SSM based on the Seebeck effect may be used to acquire a differential wave function image based on a heat-transfer gap and a temperature bias.

Additionally, in the present disclosure, an SSM image may refer to an image acquired by the above-described SSM. Hereinafter, an operation of simulating an atomic-resolution image of a material surface corresponding to an SSM image, using a computer will be further described.

Hereinafter, a simulation of an SSM image will be described based on graphene, however, a material is not limited to graphene. Accordingly, a simulation method according to embodiments may be applicable to other materials having a predetermined atomic structure.

In the present disclosure, a Seebeck coefficient may include, for example, a position-dependent Seebeck coefficient or a coherent Seebeck coefficient S_(coh)(r), and a local Seebeck coefficient S(r) or S(r; r′).

Heat may be perceived to be diffusive and transported incoherently by charge carriers (for example, electrons and holes) and lattice vibrations (for example, phonons) in a material. Accordingly, it may be very difficult to perform local imaging of a material and electronic states of the material through heat transport. Recently, however, Cho et al. reported that a series of atomic wave function images are obtained by performing local thermoelectric imaging on a surface of epitaxial graphene using a heat-based scanning probe microscope. The above report may generate a question on how to measure an atomic variation in a unit cell in a heat transport experiment. To answer the question, fundamental physics of thermoelectricity or the Seebeck effect may need to be re-evaluated from typical length scales to atomic length scales, in addition to elucidating an imaging mechanism of a scanning thermoelectric microscope.

The present disclosure may provide a theory of a scanning thermoelectric microscope with an atomic resolution based on mesoscopic electrons and heat transport characteristics. The theory, beginning with a macroscopic general transport equation and an electrostatic equation, may show feasibility and mechanisms of an atomic scale imaging technology with thermoelectric measurements. Computer simulations of a thermoelectric image may be combined with experimental results, and may be effectively used to identify atomic scale defects in graphene.

FIG. 1 schematically illustrates an atomic resolution of a scanning thermoelectric microscope according to an embodiment.

As shown in FIG. 1, a probe 10 of a scanning thermoelectric microscope at a first temperature T₁ may be in contact with graphene 20 at a second temperature T₂. Measurement may be performed through a modified ultrahigh vacuum contact mode atomic force microscope (AFM). For convenience, in the present disclosure, graphene is described, however, a simulation method according to an embodiment of the present invention may not be interpreted to be limited to graphene. A temperature gradient ∇T may be formed in the vicinity of a tip-sample contact area by a difference between the first temperature T₁ and the second temperature T₂, and an effective temperature drop ΔT_(coh)(r) may be induced in an interface between a tip and a sample. Accordingly, a thermoelectric voltage showing a phase difference of 180° with an atomic corrugation z(r) may be accurately measured by a high-impedance voltmeter in an actual scanning thermoelectric microscope. A measured local thermoelectric voltage may be expressed as shown in Equation 1 below. V(r)=V _(diff) +S _(coh)(r)ΔT _(coh)(r)  [Equation 1]

In Equation 1, V_(diff) denotes a thermoelectric voltage drop in a diffusive transport region in both the tip and the sample, S_(coh)(r) and ΔT_(coh)(r) denote a position-dependent Seebeck coefficient and an effective temperature drop at the interface between the tip and sample, respectively. Electrons and heat may be coherently transported. A coherent thermoelectric voltage V(r), represented by Equation 1, may correspond to an atomic resolution observed in a scanning thermoelectric microscope, which will be described below.

When the temperature gradient ∇T exists in a macroscopic electro-conductive system, electrons or charged particles may be transported based on an electrostatic field E and a driving force for particle diffusion under the temperature gradient ∇T. An electric current density J(r) at a local site may be expressed by a general transport equation J(r)=σ [E(r)−S(r)∇T(r)] in which σ denotes an electrical conductivity, and S(r) denotes a local Seebeck coefficient or thermopower. For example, an ideal voltmeter may satisfy an open-circuit limit, that is, J(r)=0 and accordingly, charged particles may experience a balance between an electrostatic force and a thermopower force. A built-in potential may be expressed as shown in Equation 2 below. E(r)=S(r)∇T(r)  [Equation 2]

A temperature profile T(r) may be primarily determined by thermal transport properties of a system, for example, thermal conductivities of constituent materials and interfacial thermal conductivity between the materials. Equation 2 may be derived for macroscopic diffusive systems, however, there is no limitation thereto. Accordingly, Equation 2 may also be applied to microscopic systems. When an external electric field is not applied, E(r) may indicate only a built-in electric field resulting from a thermal diffusion-induced charge distribution q^(th)(r). A Gauss's law may be applied to the built-in electric field and a charge density, that is, ∇·E(r)=q^(th)(r). Accordingly, Equation 2 may lead to Equation 3 shown below. ∇·(S(r)∇T(r))=q ^(th)(r)  [Equation 3]

From Equation 3, the thermal diffusion-induced charge distribution q^(th)(r) may be accurately traced back from information of the local Seebeck coefficient S(r) and temperature profile T(r).

When an AFM tip is assumed to be a point probe as shown in FIG. 1, a local thermoelectric voltage may be equivalent to the “Hartree-type” electrostatic potential Γ(r)=∫q^(th)(r; r′)/|r′−r|d³r′ in which r denotes a position of a probe and V(r) denotes integration over r′ for the entire volume. Based on Equation 3, and ∇T having a value of “0” at infinity, the local thermoelectric voltage may be expressed as shown in Equation 4 below.

$\begin{matrix} {{V(r)} = {\int{{S\left( {r;r^{\prime}} \right)}{{\nabla{T\left( {r;r^{\prime}} \right)}} \cdot \frac{r^{\prime} - r}{{{r^{\prime} - r}}^{3}}}d^{3}r^{\prime}}}} & \left\lbrack {{Equation}\mspace{14mu} 4} \right\rbrack \end{matrix}$

The local thermoelectric voltage V(r) may be expressed in the form of a volume integral, not a line integral, of the local Seebeck coefficient S(r) convoluted by a radially weighted temperature gradient ∇T(r) by a factor of 1/r² in which r denotes a distance measured from a point voltage probe. Because Equation 4 is an exact expression that does not involve an approximation, Equation 4 may be generally applicable for other thermoelectric systems including non-contact STM setups.

The volume integral in Equation 4 may be split into a diffusive transport region and a coherent transport region. In the diffusive transport region, a Seebeck coefficient and a temperature profile may be determined only by material properties, for example, an electrical conductivity and a thermal conductivity. From a Mott formula, a diffusive Seebeck coefficient may be represented by

$S = {{- \frac{\left( {k_{B}\pi} \right)^{2}T}{3e}}\frac{\partial}{\partial E}\left( {\ln\left\lbrack {\sigma(E)} \right\rbrack} \right)_{E_{F}}}$ in which k_(B) denotes a Boltzmann constant, T denotes an absolute temperature, e denotes an electron charge, σ(E) denotes an energy-dependent electrical conductivity, and E_(F) denotes a Fermi energy. Because a temperature profile T(r) may mostly vary slowly, under control of a phonon mean-free-path, an almost constant thermoelectric voltage V_(diff) may be generated in the diffusive transport region.

In the coherent transport region at a tip-sample interface, transport of electrons and phonons across a junction may be accounted for by each transmission probability and electrical and thermal conductance quanta. A coherent Seebeck coefficient S_(coh)(r) may be dependent on a tip position r that is independent of internal coordinates r′, and may be expressed based on an electron transmission probability τ(E,r) of a tip-sample junction as

${S_{coh}(r)} = {{- \frac{\left( {k_{B}\pi} \right)^{2}T}{3e}}\frac{\partial}{\partial E}\left( {\ln\left\lbrack {\tau\left( {E,r} \right)} \right\rbrack} \right)_{E_{F}}}$ from a Landauer formula.

A local coherent thermoelectric voltage may be expressed as shown in Equation 5 below.

$\begin{matrix} \begin{matrix} {{V_{coh}(r)} = {{S_{coh}(r)}{\int{{{\nabla{T\left( {r;r^{\prime}} \right)}} \cdot \frac{r^{\prime} - r}{{{r^{\prime} - r}}^{3}}}d^{3}r^{\prime}}}}} \\ {= {{S_{coh}(r)}\Delta\;{T_{coh}(r)}}} \end{matrix} & \left\lbrack {{Equation}\mspace{14mu} 5} \right\rbrack \end{matrix}$

In Equation 5, a volume integral of a weighted temperature profile may be set to an effective temperature drop ΔT_(coh)(r) that is sensitively dependent on a local geometry and a near-probe temperature profile due to 1/r². Physically, the effective temperature drop ΔT_(coh)(r) may correspond to an interfacial temperature drop that is known to exist at a thermal boundary between two different thin-film materials. Generally, an interfacial temperature drop at thin films may be determined based on vibrational spectra and interaction strength of involved materials.

For example, the effective temperature drop ΔT_(coh)(r) may be expressed as

${\Delta\;{T_{coh}(r)}} = {\frac{\Delta\; Q\; 3h}{\pi^{2}k_{B}^{2}T}\frac{1}{\tau^{ph}(r)}}$ from the Landauer formula, based on a phonon transmission probability of the tip-sample junction. The phonon transmission probability of the tip-sample junction may be represented by τ^(ph)(r).

Based on the diffusive thermoelectric voltage V_(diff) and the coherent thermoelectric voltage V_(coh)(r), a total local thermoelectric voltage may be derived from Equation 1. Accordingly, how thermopower profiling works for a p-n junction, and how local thermoelectric imaging works seamlessly from micrometer to sub-angstrom scales may be explained.

The coherent Seebeck coefficient S_(coh)(r), for example, a position-dependent Seebeck coefficient may be expressed as shown in Equation 6 below.

$\begin{matrix} {{S_{coh}(r)} = {{- \frac{1}{eT}}\frac{\int{{\tau^{el}\left( {E,r} \right)}\left( {E - E_{F}} \right)\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}{\int{{\tau^{el}\left( {E,r} \right)}\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}}} & \left\lbrack {{Equation}\mspace{14mu} 6} \right\rbrack \end{matrix}$

In a weak-coupling limit, the electron transmission probability τ(E,r) may be approximated as τ(E,r)∝ N_(e) ^(tip)(E)N_(e) ^(sample)(E,r) in which N_(e)(E,r) denotes a local density of states (DOS).

From the Landauer formula, the coherent Seebeck coefficient S_(coh)(r) may be expressed as S_(coh)(r)=S^(tip)+S^(sample)(r). A sample Seebeck coefficient S^(sample)(r) may be expressed as shown in Equation 7 below.

$\begin{matrix} {{S^{sample}(r)} = {{- \frac{1}{eT}}\frac{\int{{N_{e}^{sample}\left( {E,r} \right)}\left( {E - E_{F}} \right)\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}{\int{{N_{e}^{sample}\left( {E,r} \right)}\left( {- \frac{\partial f}{\partial E}} \right){\mathbb{d}E}}}}} & \left\lbrack {{Equation}\mspace{14mu} 7} \right\rbrack \end{matrix}$

In Equation 7, ƒ denotes a Fermi-Dirac distribution function at an absolute temperature T. For example, at 300 kelvin (K), (−∂ƒ ∂E) may act as an integration window of ˜0.1 electronvolt (eV) near the Fermi energy E_(F). The sample Seebeck coefficient S^(sample)(r) may be either positive or negative depending on asymmetry or energy derivative of the local DOS N_(e) ^(sample)(E,r) at the Fermi energy E_(F).

Because the local DOS N_(e) ^(sample)(E,r) may be obtained from first-principles calculations for a material surface, Equation 1 may serve as a foundation of a thermoelectric image simulation when the effective temperature drop ΔT_(coh)(r) is known at the tip-sample junction. Generally, information of the effective temperature drop ΔT_(coh)(r) may not be known either from an experiment or simulation. Instead, the effective temperature drop ΔT_(coh)(r) may be deduced by inverting Equation 1 with an experimental thermoelectric voltage Γ^(expt)(r) and a theoretical Seebeck coefficient S_(coh) ^(theory)(r), for example the sample Seebeck coefficient S^(sample)(r), as shown in FIGS. 2A through 2F.

FIGS. 2A through 2F illustrate an experimental thermoelectric voltage image for a defect-free region, a Seebeck coefficient image, line profiles, and a correlation between an effective temperature drop ΔT_(coh)(r) and van der Waals (vdW) energy E^(vdw)(r) according to an embodiment.

FIG. 2A illustrates an experimental thermoelectric voltage image for a defect-free region in bilayer graphene on Silicon carbide(0001) (SiC(0001)). In FIG. 2A, a pattern originated from 6√3×6√3R30° of graphene reconstructed by a SiC substrate may be observed.

FIG. 2B illustrates an enlarged area of the experimental thermoelectric voltage image of FIG. 2A. Experimental thermoelectric voltages may be sampled in a graphene unit cell, marked with a parallelogram.

FIG. 2C illustrates a computer-simulated Seebeck coefficient image of n-doped free-standing graphene. In FIG. 2C, a center void and a carbon atom may be marked with “V” and “C,” respectively.

FIG. 2D illustrates line profiles of an experimental thermoelectric voltage, a Seebeck coefficient, and vdW energy, along a line of V-C-C-V in FIG. 2C.

FIG. 2E illustrates a correlation between a deduced effective temperature drop ΔT_(coh)(r) and the vdW energy E^(vdw)(r). Referring to FIGS. 8A through 8D, a locally averaged sample Seebeck coefficient S^(sample)(r) with a disk radius of 0.5 angstrom (Å) may be used to deduce the effective temperature drop ΔT_(coh)(r). For example, the correlation between the deduced effective temperature drop ΔT_(coh)(r) and the vdW energy E^(vdw)(r) may be represented by ΔT_(coh)(r)=ƒ(E^(vdw)(r)). A function ƒ(E^(vdw)(r)) may include a linear fitting formula or an exponential function. In FIG. 2E, a linear fitting formula ΔT_(coh)(r)=αE^(vdw)(r)+β may be ΔT_(coh)(r)=118.05+1.6454 E^(vdw)(r).

FIG. 2F illustrates a reconstructed image of a theoretical thermoelectric voltage for n-doped free-standing graphene.

FIGS. 2A and 2B illustrate a measured thermoelectric voltage Γ^(expt)(r) for a defect-free region in bilayer graphene on SiC. FIG. 2C illustrates a theoretical Seebeck coefficient S^(sample)(r) calculated from Equation 7 for the n-doped free-standing graphene. Because a diamond-coated AFM tip and a graphene sample may interact through vdW interaction, the sample Seebeck coefficient S^(sample)(r) may be evaluated at a vdW equilibrium height calculated with a Lennard-Jones potential for pristine graphene, as shown in FIGS. 5A through 5D. The experimental thermoelectric image and theoretical thermopower images in FIGS. 2B and 2C may show a reasonable degree of correspondence in a pattern. In particular, a center void of carbon hexagon may exhibit more negative signals in both experimental and theoretical results than carbon atom sites. Accordingly, an SSM is characteristically different from an STM that picks up more current at charge-abundant atomic sites. However, it may be found from line profiling in FIG. 2D that the sample Seebeck coefficient S^(sample)(r) alone may not reproduce an atomic corrugation observed in the thermoelectric voltage Γ^(expt)(r).

To identify a role of a temperature at a length scale of coherent transport, Equation 1 may be inverted, and the effective temperature drop ΔT_(coh)(r) may be deduced. An almost linear correlation between the deduced effective temperature drop ΔT_(coh)(r) and the vdW energy E^(vdw)(r) may be shown in FIG. 2E. The above linear correlation may exist because the effective temperature drop ΔT_(coh)(r) may sensitively depend on inter-atomic thermal coupling at the tip-sample interface. In an example of weak thermal coupling on a carbon atom site, a large interfacial temperature drop may be expected. In an example of strong thermal coupling on a central void of carbon hexagon, a small interfacial temperature drop may be expected.

As a result, the deduced effective temperature drop ΔT_(coh)(r) may exhibit an atomic variation at a sub-angstrom scale from a carbon atom site to a hexagonal void. Therefore, the atomic resolution in the local thermoelectric voltage may be originated from an atomic variation in coherent electron transport via a Seebeck coefficient, and an atomic variation in coherent thermal transport via a thermal conductivity. FIG. 2F illustrates a reconstructed image of a theoretical thermoelectric voltage of pristine free-standing graphene, simulated with the sample Seebeck coefficient S^(sample)(r) and the linear fitting formula of the effective temperature drop ΔT_(coh)(r).

FIGS. 3A through 3D illustrate images of a thermoelectric voltage for a point defect, a topographic image, and profiles according to an embodiment.

FIG. 3A illustrates a large-area scanning image of a thermoelectric voltage for a point defect in bilayer graphene on SiC.

FIG. 3B illustrates a small-area scanning image of a thermoelectric voltage in a rotated view for the point defect in FIG. 3A. FIG. 3C illustrates a topographic image. The small-area scanning image of the thermoelectric voltage, and the topographic image may be simultaneously obtained.

FIG. 3D illustrates height profiles along a dashed line 1 and a bold line 2 of FIG. 3C. In FIG. 3D, comparing the dashed line 1 to the bold line 2, a single carbon atom site may be clearly defective along the dashed line 1.

Based on a linear correlation between the effective temperature drop ΔT_(coh)(r) and the vdW energy E^(vdw)(r), atomic scale defects on a surface of graphene may be identified by comparing experimental thermoelectric images to simulated thermoelectric voltage images. FIGS. 3A and 3B show thermoelectric images of the point defect in bilayer graphene on SiC. A topographic analysis in FIGS. 3C and 3D may be used to infer that the point defect in an experiment is associated with a single carbon atom site.

FIGS. 4A through 4F illustrate simulated images of Seebeck coefficients and thermoelectric voltages near point defects in n-doped free-standing graphene according to an embodiment.

FIGS. 4A and 4D illustrate a simulated image of a Seebeck coefficient and a simulated image of a thermoelectric voltage for a point defect of a single carbon vacancy (V_(c)), respectively.

FIGS. 4B and 4E illustrate a simulated image of a Seebeck coefficient and a simulated image of a thermoelectric voltage for a point defect of substitutional nitrogen (N_(c)), respectively.

FIGS. 4C and 4F illustrate a simulated image of a Seebeck coefficient and a simulated image of a thermoelectric voltage for a point defect of a defect complex (V_(c)-O_(c)) of a carbon vacancy and substitutional oxygen, respectively.

In FIGS. 4A through 4F, atomic models may be shown in an inset.

To map a defect with an atomic model, a simulation may be performed on thermoelectric images of the V_(c) of FIGS. 4A and 4D, the N_(c) of FIGS. 4B and 4E, and the V_(c)-O_(c) of FIGS. 4E and 4F in the n-doped free-standing graphene.

Because an electron of the V_(c) defect is located below the Fermi energy as shown in FIG. 7A through 7D, thermopower and thermoelectric images of the V_(c) defect may have a bright region at a defect site, as opposed to an experimental image. On the other hand, because electrons of the N_(c) defect and the V_(c)-O_(c) defect are located above the Fermi energy as shown in FIGS. 8A through 9D, thermopower and thermoelectric images of the N_(c) defect and the V_(c)-O_(c) defect may have a dark region, similarly to the experimental image. While the size of the dark region matches well with a size of the N_(c), a symmetry of the image may match well with an image of the V_(c)-O_(c). Atomic oxygen may exist during sample growth and formation energy of the V_(c)-O_(c) may be sufficiently smaller than that of the V_(c) and accordingly, a point defect in an experiment may be highly likely to be V_(c)-O_(c). A simulation result including positive signals at V_(diff) of −0.6 millivolt (mV) as shown in FIGS. 10A through 10D may reproduce several conspicuous features, compared to an experimental image. Indeed, positive thermoelectric signals near a defect site in the experiment may be associated with fine details of wave function overlap, atom-by-atom thermal coupling, a substrate effect, and diffusive shift V_(diff) that may be limitedly reflected in a current simulation scheme.

Comparing a heat-based scanning thermoelectric microscopy to an STM may be important. Both the above techniques may share common features and functionalities as types of a scanning probe microscopy that provides real-space images of wave functions. The STM may measure a tunneling current by applying a voltage drop across a vacuum-tunneling gap, whereas the heat-based scanning thermoelectric microscopy may measure voltage differences by applying a temperature drop across an interface or across a heat-transfer gap. As a result, Fermi electrons may be only perturbed at the first order by a temperature bias in the heat-based scanning thermoelectric microscopy, in contrast to the zero-order perturbation by a voltage bias in the STM. Accordingly, a scanning thermoelectric microscopy may be useful to differentially analyze Fermi electronic states, even at room temperature, unlike the STM.

Hereinafter, the first-principles calculations will be further described.

A mesoscopic Seebeck coefficient in Equation 7 may be obtained using Kohn-Sham wave functions and a local DOS from first principles density-functional theory calculations of graphene. To calculate a ground-state total-energy, VASP software may be used. A “12×12” graphene supercell may be used for modeling the point defects, that is, the V_(c), N_(c), and V_(c)-O_(c) in graphene. For example, when a Dirac point is set to 0 eV, Fermi energy E_(F) of 0.3 eV may be used to calculate a local DOS. At an absolute temperature T of 315 K, the sample Seebeck coefficient S^(sample)(r) in Equation 7 may be calculated.

Hereinafter, atomically-varying effective temperature drops will be described.

To deduce an effective temperature drop ΔT_(coh)(r) from an experimental thermoelectric voltage Γ^(expt)(r) and a theoretical Seebeck coefficient S^(sample)(r) by inverting Equation 1, a locally-averaged Seebeck coefficient S^(sample)(r) within a specific disk radius R may be used. This is a method based on an effect of a finite size of a probe. In FIGS. 2E and 2F, a locally-averaged Seebeck coefficient S^(sample)(r) with a disk radius R of 0.5 Å may be used. For comparison, an original Seebeck coefficient S^(sample)(r) and a locally-averaged Seebeck coefficient S^(sample)(r) with a disk radius R of 0.3 Å may be used, as shown in FIGS. 6A through 6D.

The deduced effective temperature drop ΔT_(coh)(r) shows an almost linear correlation or an almost exponential correlation with vdW energy E^(vdw)(r) that is calculated with Lennard-Jones 12-6 potentials, as shown in Equation 8 below.

$\begin{matrix} {{E^{vdW}\left( r_{i} \right)} = {\sum\limits_{j}{4{ɛ_{ij}\left\lbrack {\left( \frac{\sigma_{ij}}{{r_{i} - r_{j}}} \right)^{12} - \left( \frac{\sigma_{ij}}{{r_{i} - r_{j}}} \right)^{6}} \right\rbrack}}}} & \left\lbrack {{Equation}\mspace{14mu} 8} \right\rbrack \end{matrix}$

In Equation 8, r_(i) and r_(j) denote an atomic position of a tip, an atomic position of a sample, respectively, and ε and σ denote Lennard-Jones parameters. The tip may be modeled with a single carbon atom. Parameters ε_(ii) and σ_(ii) for a carbon atom, a nitrogen atom, and an oxygen atom may be listed as shown in Table 1, and ε_(ij)=∇(ε_(ii)ε_(jj)) and σ_(ij)=(σ_(ii)+σ_(jj))/2 may be used. The vdW energy may be summed when an atom-atom distance is less than 15 Å, and an equilibrium height z(r) at minimum vdW energy may be used to evaluate a Seebeck coefficient S^(sample)(r).

TABLE 1 i ε_(ii) (meV) σ_(ii) (Å) C 4.20 3.37 N 7.41 3.25 O 9.12 2.96

Hereinafter, a statistically defined Fermi temperature will be further described.

In-depth physics of the sample Seebeck coefficient S^(sample)(r) expressed by Equation 7 may be appreciated. The sample Seebeck coefficient S^(sample)(r) may be expressed as shown in Equation 9 below.

$\begin{matrix} {{S^{sample}(r)} = {{- \frac{\left( {k_{B}\pi} \right)^{2}T}{3e}}\frac{\partial}{\partial E}\left( {\ln\left\lbrack {N_{e}^{sample}\left( {E,r} \right)} \right\rbrack} \right)_{E_{F}}}} & \left\lbrack {{Equation}\mspace{14mu} 9} \right\rbrack \end{matrix}$

When k_(B) ln [N_(e) ^(sample)(E,r)] is defined as an electronic entropy, energy derivative of an entropy may correspond to a reciprocal temperature, expressed as shown in Equation 10 below.

$\begin{matrix} {\frac{1}{T_{F}(r)} = {\frac{\partial}{\partial E}\left( {k_{B}{\ln\left\lbrack {N_{e}^{sample}\left( {E,r} \right)} \right\rbrack}} \right)_{E_{F}}}} & \left\lbrack {{Equation}\mspace{14mu} 10} \right\rbrack \end{matrix}$

In Equation 10, T_(F)(r) denotes a position-dependent statistically-defined Fermi temperature. A newly defined Fermi temperature T_(F)(r) of a real material may be a generalized version of a “Fermi temperature” T_(F)=E_(F)/k_(B) that is only applicable to a three-dimensional (3D) free electron gas model. Additionally, a Seebeck coefficient may be expressed as shown in Equation 11 below.

$\begin{matrix} {{S(r)} = {{- \;\frac{\pi^{2}}{3}}\frac{k_{B}}{e}\frac{T}{T_{F}(r)}}} & \left\lbrack {{Equation}\mspace{14mu} 11} \right\rbrack \end{matrix}$

Equation 11 may be interpreted as a ratio of a thermal equilibrium temperature T and Fermi temperature T_(F) of a material. A statistically-defined Fermi temperature may be a material property, not a real temperature, and may be either positive or negative depending on a slope of N_(e) ^(sample)(E,r). The Seebeck coefficient in Equation 11 may be conceptually connected to other electron-related thermal characteristics of a material, for example, an electronic heat capacity c_(v)=π²N_(e)(E_(F))k_(B) ²T/3 and a quantum of thermal conductance g₀=π²k_(B) ²T/3h. For example, the Seebeck coefficient and electronic heat capacity may be expressed as shown in Equation 12 below.

$\begin{matrix} {{S(r)} = {\frac{c_{v}}{e}\frac{\partial}{\partial E}\left( \frac{1}{N_{e}^{sample}\left( {E,r} \right)} \right)_{E_{F}}}} & \left\lbrack {{Equation}\mspace{14mu} 12} \right\rbrack \end{matrix}$

FIGS. 5A through 5D illustrate a model and a DOS of pristine free-standing graphene, vdW energy, and a vdW topography according to an embodiment.

FIG. 5A illustrates a ball-and-stick model of the pristine free-standing graphene.

FIG. 5B illustrates the DOS of the pristine free-standing graphene. In FIG. 5B, zero energy may indicate a charge-neutrality point or a Dirac point that is denoted by a bold dashed line. A fine dashed line may mark a Fermi energy of 0.3 eV used for a thermoelectric simulation of FIGS. 4A through 4F.

FIG. 5C illustrates a computer-simulated image of vdW energy E^(vdw)(r).

FIG. 5D illustrates a computer-simulated image of vdW topography z(r) at minimum energy.

FIGS. 6A through 6F illustrate locally-averaged Seebeck coefficients, a correlation between vdW energy and an effective temperature drop, and thermoelectric voltages according to an embodiment.

FIG. 6A illustrates a computer-simulated image of a locally-averaged Seebeck coefficient within a disk radius R of 0 Å.

FIG. 6B illustrates a computer-simulated image of a locally-averaged Seebeck coefficient within a disk radius R of 0.3 Å.

FIG. 6C illustrates a correlation between vdW energy and an effective temperature drop deduced with the disk radius R of 0 Å.

FIG. 6D illustrates a correlation between vdW energy and an effective temperature drop deduced with the disk radius R of 0.3 Å.

FIG. 6E illustrates a reconstructed image of a thermoelectric voltage for the disk radius R of 0 Å.

FIG. 6F illustrates a reconstructed image of a thermoelectric voltage for the disk radius R of 0.3 Å.

FIGS. 7A through 7D illustrate a model and a DOS of defective free-standing graphene, vdW energy, and a vdW topography according to an embodiment.

FIG. 7A illustrates a ball-and-stick model of defective free-standing graphene with a single V_(c).

FIG. 7B illustrates a DOS of defective free-standing graphene with a single V_(c), that is, a local DOS of dangling-bonded carbon C₃. In FIG. 7B, zero energy may indicate an original Dirac point, a bold dashed line may indicate original Fermi energy of defective graphene, and a fine dashed line may indicate the elevated Fermi energy of 0.3 eV used for a thermoelectric simulation of FIGS. 4A and 4D.

FIG. 7C illustrates a computer-simulated image of vdW energy E^(vdw)(r).

FIG. 7D illustrates a computer-simulated image of a vdW topography z(r) at minimum energy.

FIGS. 8A through 8D illustrate a model and a DOS of defective free-standing graphene with a single N_(c), vdW energy, and a vdW topography according to an embodiment.

FIG. 8A illustrates a ball-and-stick model of defective free-standing graphene with a single N_(c).

FIG. 8B illustrates a DOS of defective free-standing graphene with a single N_(c), marked with a medium-thick line. In FIG. 8B, a local DOS for a nitrogen atom and a local DOS of three neighboring carbon atoms may be displayed. In FIG. 8B, Zero energy may indicate an original Dirac point, a bold dashed line may indicate original Fermi energy of defective graphene, and a fine dashed line may indicate Fermi energy of 0.3 eV used for a thermoelectric simulation of FIGS. 4B and 4E.

FIG. 8C illustrates a computer-simulated image of vdW energy E^(vdw)(r).

FIG. 8D illustrates a computer-simulated image of a vdW topography z(r) at minimum energy.

FIGS. 9A through 9D illustrate a model and a DOS of defective free-standing graphene with a V_(c)-O_(c), vdW energy, and a vdW topography according to an embodiment.

FIG. 9A illustrates a ball-and-stick model of the defective free-standing graphene with the V_(c)-O_(c). In FIG. 9A, an oxygen atom may be marked with bold lines.

FIG. 9B illustrates a DOS of the defective free-standing graphene with the V_(c)-O_(c). FIG. 9B may show a local DOS for atoms around a vacancy. In FIG. 9B, zero energy may indicate an original Dirac point, a bold dashed line may indicate original Fermi energy of defective graphene, and a fine dashed line may indicate Fermi energy of 0.3 eV used for a thermoelectric simulation of FIGS. 4C and 4F.

FIG. 9C illustrates a computer-simulated image of vdW energy E^(vdw)(r).

FIG. 9D illustrates a computer-simulated image of a vdW topography z(r) at minimum energy.

FIGS. 10A through 10D illustrate a simulated thermoelectric voltage image for a V_(c)-O_(c), a fast Fourier transform (FFT) image corresponding to the simulated thermoelectric voltage image, an experimental thermoelectric image, and an FFT image corresponding to the experimental thermoelectric image according to an embodiment.

FIG. 10A illustrates a computer-simulated image of a thermoelectric voltage for the V_(c)-O_(c) at V_(diff) of −0.6 mV.

FIG. 10B illustrates an FFT image corresponding to the computer-simulated image of FIG. 10A. In FIG. 10B, an arrow A may correspond to a reciprocal lattice, and an arrow B may be caused by intervalley scattering of a Fermi wave vector k_(F).

FIG. 10C illustrates an experimental thermoelectric image that may correspond to the large-area scanning image of the thermoelectric voltage of FIG. 3A.

FIG. 10D illustrates an FFT image corresponding to the experimental thermoelectric image of FIG. 10C. In FIG. 10D, an arrow A may correspond to a reciprocal lattice, and an arrow B may be caused by intervalley scattering of a Fermi wave vector k_(F). Additionally, an arrow C may represent a pattern originated from 6√3×6√3R30° reconstruction of a surface of SiC(0001).

Although a few exemplary embodiments of the present invention have been shown and described, the present invention is not limited to the described exemplary embodiments. Instead, it would be appreciated by those skilled in the art that changes may be made to these exemplary embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents. 

What is claimed is:
 1. A computer-aided simulation method for a scanning Seebeck microscope (SSM) image, wherein a computer calculates a local thermoelectric voltage for a position of a voltage probe in a material surface, to acquire the SSM image corresponding to the position, using the following equation: ${V(r)} = {\int{{S\left( {r;r^{\prime}} \right)}{{\nabla{T\left( {r;r^{\prime}} \right)}} \cdot \frac{r^{\prime} - r}{{{r^{\prime} - r}}^{3}}}d^{3}r^{\prime}}}$ in which V(r) denotes the local thermoelectric voltage, S(r; r′) denotes a local Seebeck coefficient, ∇T(r; r′) denotes a temperature gradient radially weighted by a factor of 1 r², r denotes a position vector measured from a point voltage probe, r′ denotes material internal coordinates, |r′−r| denotes an absolute value of r′−r, and d³r′ denotes a volume integral.
 2. A computer-aided simulation method for a scanning Seebeck microscope (SSM) image, wherein a computer calculates a thermoelectric voltage based on a Seebeck coefficient and a change in a temperature between a probe and a surface of a material, the thermoelectric voltage being used to display an image of a predetermined local region of the material.
 3. The computer-aided simulation method of claim 2, wherein the Seebeck coefficient is calculated by the computer based on at least one of an electron charge, an absolute temperature, a Fermi energy, and a Fermi-Dirac distribution function.
 4. The computer-aided simulation method of claim 2, wherein the Seebeck coefficient is calculated by the computer based on at least one of an electrical conductivity, an electron transmission probability, and a local density of states (DOS).
 5. The computer-aided simulation method of claim 2, wherein the Seebeck coefficient comprises a unique characteristic of the material.
 6. The computer-aided simulation method of claim 2, wherein the change is calculated by the computer based on at least one of a thermal conductivity, a phonon transmission probability, and an interaction between the probe and the material.
 7. The computer-aided simulation method of claim 2, wherein the material comprises graphene.
 8. The computer-aided simulation method of claim 2, wherein the thermoelectric voltage is induced by a temperature difference between a first temperature of the probe and a second temperature of the material. 